Equilibrium and non-equilibrium dynamics of the sub-ohmic spin-boson model 
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Employing the non-perturbative numerical renormalization group method, we study the dynamics 
of the spin-boson model, which describes a two-level system coupled to a bosonic bath with spectral 
density J {to) oc u) s . We show that, in contrast to the case of ohmic damping, the delocalized phase 
of the sub-ohmic model cannot be characterized by a single energy scale only, due to the presence of 
a non-trivial quantum phase transition. In the strongly sub-ohmic regime, s <C 1, weakly damped 
coherent oscillations on short time scales are possible even in the localized phase - this is of crucial 
relevance, e.g., for qubits subject to electromagnetic noise. 
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Models of quantum dissipation 0, 0] have gained sig- 
nificant attention over the last years, due to their wide 
range of applications from the effect of friction on the 
electron transfer in biomolecules Q to the description of 
the quantum entanglement between a qubit and its envi- 
ronment [J, [5( (for further applications see Refs. yj, |2|) . 

The simplest quantum-dissipative models belong to the 
class of impurity models and consist of a small quantum 
system coupled to a bath of harmonic oscillators. One 
familiar bosonic impurity model is the spin-boson model, 
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It describes a generic two-level system, represented by 
the Pauli matrices Oj , which is linearly coupled to a bath 
of harmonic oscillators, with creation operators a\ and 
frequencies cjj. The bare tunneling amplitude between 
the two spin states | j), | J.) is given by A, and e is an 
additional bias (which is zero in the following, except for 
the preparation of the initial state as discussed below). 
The coupling between spin and bosonic bath is specified 
by the bath spectral function J{uj) = TtJ2i ^H^i^ ~ ^i)- 
The asymptotic low-temperature behavior is determined 
by the low-energy part of the spectrum. Discarding high- 
energy details, the standard parametrization is 

J(uj) — 2irau 1 c ~ s uj s , < u) < uj c , s > -1 (2) 

where the dimensionless parameter a characterizes the 
dissipation strength, and u> c is a cutoff energy. 

In case of ohmic dissipation, s = 1, a quantum tran- 
sition of Kosterlitz-Thouless type separates a localized 
phase at a > a c , displaying a doubly degenerate ground 
state, from a delocalized phase at a < a c with a unique 
ground state [l], Q • The delocalized regime is character- 
ized by a finite effective tunnel splitting, A r , between the 
two levels, whereas the tunnel splitting renormalizes to 
zero in the localized phase. For A <IC lu c the transition 
occurs at a c = 1. 

The sub-ohmic case, s < 1, d, [ESU, turns out to 
be different. For A/oj c — > the system is localized for 
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FIG. 1: Zero-tem per ature phase diagram of the spin-boson 
model from NRG Q, for fixed A = 0.1 andw c = 1. The 
axes denote the bath exponent s and the dissipation strength 
a. The circles indicate the parameter values for which results 
will be presented in this paper. 



any non-zero coupling 0, Q , but for large A a delocal- 
ized phase was argued to exist [1, 0] ■ We have recently 
8] shown, using an extension of the non-perturbative 
Numerical Renormalization Group (NRG), that a con- 
tinuous quantum phase transition (QPT) occurs for all 
< s < 1, in contrast to earlier proposals Q. The nu- 
merically determined equilibrium phase diagram is shown 
in Fig. [B 

Detailed studies of non-equilibrium properties, like the 
quantity P(t) = (a z (t)} for a spin initially prepared in the 
|) state, have mainly covered the ohmic damping case 
for small A [l], [l(| [ll[ . Here, weakly damped oscillations 
can only be observed for a < 1/2, whereas 1/2 < a < 1 
leads to overdamped behavior, i.e., exponential decay of 
P(t) to zero. Finally, in the localized phase, a > 1, P(t) 
decays to a finite value. Studies of the dynamics in the 
sub-ohmic case have so far been restricted to the use of 
perturbative methods [ijj]- Considering that these miss 
the QPT and the localized phase present for s < 1, the 
validity of the corresponding results is questionable. 

The purpose of this paper is to provide non- 
perturbative results for the dynamics of the sub-ohmic 
spin-boson model, for the whole range of model param- 
eters. Our main results can be summarized as follows: 
Whereas the delocalized phase of the ohmic model is 
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dominated by a single energy scale, A r , only the sub- 
ohmic model is more complicated due to the presence 
of a second-order QPT. Near the transition a quantum 
critical (QC) crossover scale T* appears 0], which - in 
the strongly sub-ohmic regime - co-exists with a larger 
scale A r that can still be identified with the renormal- 
ized tunnel splitting. As we demonstrate below, T* van- 
ishes at the transition to the localized phase, whereas 
A r stays finite. The latter fact is crucial for the non- 
equilibrium dynamics: we find the presence of coherent 
weakly damped oscillations, with frequency A r , even in 
the localized phase, a > a c , for s < 1. This is partic- 
ularly relevant for the case s = 0, related to so-called 
1/f noise in an electromagnetic environment. Here, the 
equilibrium spin-boson model is always in the localized 
phase, however, coherent oscillations are still possible for 
small a and short times. 

NRG. We employ Wilson's NRG method [13] to study 
the sub-ohmic spin-boson model, utilizing two recently 
developed extensions: (i) Refs.l8l.fl4l generalized the NRG 
at thermodynamic equilibrium to impurity models with 
a bosonic bath, (ii) Ref. [TBI proposed an algorithm to 
study non-equilibrium dynamics in real time and applied 
it to the Anderson and Kondo models. (For an earlier 
non-equilibrium NRG approach see Ref. LL6|.) The time- 
dependent NRG provides a spectral representation of the 
time-independent Hamiltonian , governing the time 
evolution for t > 0, at all energy scales using a complete 
basis set, and expresses the real-time dynamics of ob- 
servables in terms of a summation over reduced density 
matrices [IB] ] . The latter contain all information on deco- 
herence and dissipation. Two independent bosonic NRG 
runs [14j are required, one for the initial density matrix 
and the other for the approximate eigenbasis of H f . In 
order to accurately simulate the continuum limit with a 
NRG chain of finite length, we average over Nz different 
bath discretizations for a fixed discretization parameter 
A. All numerical results below are for temperature T = 0, 
unless otherwise noted. 

Renormalization group flow and crossovers. To set the 
stage, we summarize the renormalization group flow of 
the spin-boson model [1, [3, [l7|. In the following, the 
terms localized (delocalized) are defined through the im- 
purity entropy [3, EH being In 2 (zero) . Note that this 
has to be contrasted with localization in the sense that 
a system initially prepared with the impurity spin in one 
specified direction remains in this spin state under time 
evolution. For any finite temperature, thermal excita- 
tions destroy localization in this sense (see Ref. El). For 
a discussion of this point, in particular the connection 
between NRG flow and thermodynamic properties, see 
Sees. Ill and IV in Ref. El- 

In the ohmic case, the flow for a < a c is from the 
localized towards the delocalized phase (upon lowering 
the energy), with A r being the crossover energy scale. 
Thus the behavior for energies or temperatures below 
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FIG. 2: Schematic crossover diagrams for the spin-boson 
model for different bath exponents s, deduced from NRG, 
with regimes of localized, delocalized, and quantum critical 
(QC) behavior. The open dots denote QPTs, for details see 
text. (Note that the a scale is different for the three panels.) 



A r is delocalized, Fig. [2k- 

For < s < 1 the critical fixed point controls the 
physics of the QC region. It is bounded by crossover 
lines T* cx \a — a c \" and covers a large portion of the 
phase diagram due to the large value of the correlation 
length exponent v, Fig. [2}d. (In fact, v diverges both as 
s — > I ~~ and s — > + [l7|.) The critical fixed point merges 
with the localized (delocalized) one as s — > l - (s — > + ), 
implying that the characteristics of the QC regime be- 
comes more "delocalized" as s is decreased towards zero. 
For s not too close to I there also is a distinct crossover 
from the high-temperature localized regime (T 3> A) to 
the delocalized or critical regimes. 

Consequently, for s = and small a the flow is first 
from localized to delocalized, with a crossover scale A r , 
finally to the localized fixed point (which controls the 
ground state for any a) - the latter crossover is charac- 
terized by the scale T*, with T* ~ uj c exp(— A/au> c ). As 
shown in Fig. [2J;, the low-energy crossover at s = is 
thus opposite to the one at s = 1! 

In general, we may expect coherent, weakly damped 
dynamics, with a rate A r , in the delocalized regime (and 
in the QC one for small s as well). 

Equilibrium dynamics. Let us focus on the Fourier 
transform C(uj) of the symmetrized equilibrium correla- 
tion function C(t) = i([er 2 (i), cr z } + ). We start by noting 
the low-frequency asymptotics: In the delocalized phase 
C{uj) cx oj s @, [3], whereas in the localized phase C{lo) 
displays a A5{lo) contribution, which reflects the doubly 
degenerate ground state, and A plays the role of an order 
parameter (equivalent to the thermodynamic expectation 
value (o~ z )). Turning to the finite-frequency behavior, our 
numerics for s = 1 reproduces the well-known results (see 



Sec. V in Ref.ll4T), i.e., C{u>) is dominated by a single peak 
at A r which shifts to lower frequencies with increasing a 
and disappears at the transition, a = a c . (Results for s 
down to 0.8 are qualitatively similar.) 

For smaller s the influence of the critical fixed point 
becomes increasingly visible, see Fig. |3k for s = 0.5. 
For small a we observe that C{lo) is dominated by a 
peak close to A, as expected. However, with increasing 
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FIG. 3: Equilibrium spin correlation function C(ui) for a) 
s = 0.5 and a = 0.02,0.06,0.09,0.10,0.107,0.11; b) s = 0.1 
and a = 0.002,0.004,0.006,0.007,0.01; the thick curves are 
very close to the critical point (for s = 0.5 we find a c — 0.1065 
and s = 0.1, a c — 0.0071). The curves correspond to the 
parameter values shown as circles in Fig. [TJ (The peak width 
at small a is an artifact of the NRG broadening.) 



a weight is transferred to smaller frequencies, leaving a 
shoulder feature close to A intact. For a < a c , C(ui) 
has a pronounced peak at ui ~ T* - this peak separates 
the QC divergence C(w) oc uj~ s at intermediate energies 
from the low-energy uu s behavior. At criticality, T* — > 0, 
and C(uj) oc lu~ s down to lowest energies. 

A further decrease of the bath exponent s - see Fig. [3Jd 
for s = 0.1 - shows that the high- frequency peak in the 
vicinity of A is suppressed more slowly upon increasing 
a, reflecting the fact that the QC behavior is more "de- 
localized" for small s. In particular, this peak, which is 
usually taken as indication for the presence of coherent 
weakly damped dynamics, survives even for a > a c (!). 

Fig. [4] for s — illustrates this fact: here the ground 
state is localized for any s, but a well-defined peak at 
A r is visible for all a < 0.05 (for A/uj c = 0.1). Further, 
there is a low-energy peak at T* for small a, in addition 
to the A8(lo) contribution (not plotted) characteristic of 
the localized phase. 

Parenthetically, we note that the low-energy part of 
C(oj,T) near a = a c is expected to show universal scal- 
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FIG. 4: Equilibrium spin correlation function C(ui) for s = 
and a = 0.003, 0.005, 0.007, 0.01, 0.02. Both panels show the 
same data, but the right one has a linear u scale. The AS(lu) 
contribution is not shown. 



ing behavior, including ui/T scaling, as the critical fixed 
point obeys hyperscaling properties for all < s < 1 [17 1. 

Non-equilibrium dynamics. We now present results for 
P(£), obtained with the non-equilibrium generalization 
of the NRG [HI]. The spin is prepared in the ||) state, 
by calculating the equilibrium density matrix of the ini- 
tial Hamiltonian H % with a large spin polarization energy 
e/uj c = 100 and A = 0. At t = 0, e is switched off, and 
A is set to A/lu c = 0.1. This defines H$ , governing the 
time evolution of the initial density operator. 

Fig. [5]shows P(t) for s — 0.5, A/cj c = 0.1, and various 
values of a. For weak coupling, i.e. a <C a c , a damped 
oscillatory behavior is found. Increasing a suppresses 
the oscillations at longer time scales but maintains the 
initial ones. A very shallow oscillation at short times can 
be seen even for a > a c which is consistent with the C(u) 
data presented in Fig. [3l For a > a c P(t) shows only a 
very weak time dependence for intermediate times and is 
expected to approach a finite value P(oo) for t — > oo, an 
indication of localization (see below). 

Let us now turn to s = 0, Fig. [6] Even though any fi- 
nite value of a places the model in the localized phase in 
equilibrium, the non-equilibrium dynamics of P(t) clearly 
exhibits oscillatory behavior for small a. Roughly speak- 
ing, the time evolution mimics the RG flow and thus 
corresponds to a reduction of temperature in the phase 
diagram Fig. [2t: At short times the physics is governed 
by the delocalized fixed point; it crosses over to the lo- 
calized fixed point at long time scales. Upon increas- 
ing a the time scale of the crossover from oscillatory to 
damped behavior in P(t) is reduced, consistent with the 
low-temperature crossover line in Fig. [2J;. 

For small damping, the decay of the P(t) oscillations is 
exponential. Comparing the decay rate with the popular 
Bloch-Redfield approach [20( gives deviations of less than 
10% for s — 0.1, a — 0.002, but the agreement becomes 
worse for larger a or smaller s. For all s < 1 we observe 
that with increasing a the frequency of the initial oscil- 
lations in P(t) first decreases (consistent with the result 
of straightforward perturbation theory) , but for larger a 
and small s then slightly increases with a, indicating that 
one leaves the perturbatively accessible regime. [C(u>) in 
Figs. [3l [4] is inconclusive w.r.t. the peak position due to 
the logarithmic broadening with the NRG.] 

The long-time limit P(oo) measures the "degree of lo- 
calization", equivalent to the order parameter (a z ) in 
equilibrium [and proportional to the ^-coefficient of the 
S(u) contribution to C(uu)\. P(oo) vanishes continuously 
at the transition for any < s < 1. [For the smallest 
value of a shown in Fig. [6] for s = (a — 0.003), the 
P(t) oscillates around P(oo) w 0.1.] Note the difference 
to the ohmic case, s = 1, where P(oo) jumps from a finite 
value for a > a c to zero a < a c . 

A few remarks are in order. On general grounds, one 
may expect a power-law decay in the long-time limit of 
P(t) at the critical point, a — a c , for all < s < 1. This 
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FIG. 5: Time evolution of P(t) for s = 0.5, A/lu c 
various values of a (a c = 0.1065). 



0.1, and 



power law will be cutoff away from criticality, on a time 
scale <xl/T*. At present, the accuracy of the numerical 
data in non-equilibrium is not sufficient to verify this. 
Since the short-time behavior of P(t) is dominated by 
the high-energy properties, i.e., by A r , small but finite 
temperatures do not alter the response on time scales 
shown in Figs. [5] and [6] At finite T we only observe 
a decay of P(t) to zero for long times and parameters 
corresponding to the localized phase (not shown). 

Summary. We have studied the dynamics of the sub- 
ohmic spin-boson model, both in equilibrium and non- 
equilibrium, using non-perturbative numerical renormal- 
ization group techniques. The model displays a continu- 
ous QPT for all bath exponents < s < 1, and this leads 
to highly non-trivial dynamical properties. In contrast to 
the ohmic situation, the sub-ohmic case cannot be char- 
acterized by a single energy scale only. This is particu- 
larly striking for s<l: while the low-energy, long-time 
behavior is dominated by the presence of nearly critical 
fluctuations over a large regime of parameters, the behav- 
ior at short times or elevated energies derives from weakly 
renormalized tunneling, which can be coherent even in 
situations with a localized thermodynamic ground state. 

Our results imply that perturbative methods, only cap- 




turing the renormalized coherent tunneling, may be ap- 
plied to the sub-ohmic model at short times, but clearly 
fail in the long-time limit. Furthermore, approximations 
aiming on a description of the spin-boson physics in terms 
of a single energy scale A r only, as done in Refs. 0, are 
not applicable in the strongly sub-ohmic regime. Last 
not least, let us emphasize that the sub-ohmic situations 
of s = 1/2 and s = are of immediate experimental rele- 
vance, e.g. for the description of electromagnetic trans- 
mission lines [2l| and general 1// noise, respectively. 
Furthermore, s = 1/2 baths may be realized in the con- 
text of effective impurities in ultracold gases [22| ■ 
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FIG. 6: Time evolution of P(t) for s = 0, A/lu c = 0.1, and 
various values of a (a c = here). 



[1] 
[2] 

[3] 
[4] 



[5 
[6 

[7: 

[8 
[9 

[10 

[11 
[12 



[13 
[14 

[is; 

[16 
[17 

[is; 

[19 
[20 

[21 



[22 
[23 



A. J. Leggett et ai, Rev. Mod. Phys. 59, 1 (1987). 

U. Weiss, Quantum dissipative systems, World Scientific, 

Singapore (1999). 

A. Garg et al, J. Chem. Phys. 83, 4491 (1985). 

M. Thorwart and P. Hanggi, Phys. Rev. A 65, 012309 

(2001); M. J. Storcz and F. K. Wilhelm, ibid. 67, 042319 

(2003). 

D. V. Khveshchenko, Phys. Rev. B 69, 153311 (2004). 
S. Kehrein and A. Mielke, Phys. Lett. A 219, 313 (1996). 
H. Spohn and R. Diimcke, J. Stat. Phys. 41, 389 (1985). 
R. Bulla et al, Phys. Rev. Lett. 91, 170601 (2003). 
A. Chin and M. Turlakov, Phys. Rev. B 73, 075311 
(2006). 

T. A. Costi and C. Kieffer, Phys. Rev. Lett. 76, 1683 
(1996). 

M. Keil and H. Schoeller, Phys. Rev. B 63, 180302(R) 
(2001). 

A. Shnirman and G. Schon, in: Quantum Noise in 
Mesoscopic Physics, ed. Yu. Nazarov, Kluwer, Dordrecht 
(2003). 

K. G. Wilson, Rev. Mod. Ph ys. 47, 773 (1975); R . Bulla, 
T. Costi, and Th. Pruschke, |cond-mat/0701105| 
R. Bulla et ai, Phys. Rev. B 71, 045122 (2005). 

F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 
(2005) and Phys. Rev. B 74, 245113 (2006). 

T. A. Costi, Phys. Rev. B 55, 3003 (1997). 
M. Vojta et al, Phys. Rev. Lett. 94, 070604 (2005). 
Impurity entropy is defined as the impurity contribution 
to the thermodynamic entropy of the model |T} ■ 
H. Spohn and W. Zwerger, J. Stat. Phys. 94, 1037 (1999). 
V. May and O. Kiihn, Charge and Energy Transfer Dy- 
namics in Molecular Systems, Wiley- VCH, Berlin (2000). 

G. L. Ingold and Yu. V. Nazarov, in: Single Charge Tun- 
neling, eds. H. Grabert and M. H. Devoret, Plenum Press, 
New York (1992). 

W. Zwerger, private communication, see also: A. Re- 
cati et al., cond-mat/0212413 and Phys. Rev. Lett. 94, 
040404 (20057 - 

The NRG parameters are: A = 2, JV S = 100, N b = 8, Nw = 



5 



100. The non-equilibrium calculations employ Nz = 32 
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